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ABSTRACT 


Expressions have been obtained for the higher order Thermo elastic 
properties of the Gruneisen parameter in the limit of infinite pressure The 
usefulness of these expressions has been demonstrated by considering a 
relationship between reciprocal Gruneisen parameter and the pressure 
bulk modulus ratio at finite pressure. The results have been obtained using 
the shanker reciprocal gamma relationship. The results have also been 
determined using the Stacey-Davis formulations for the third order 
Gruneisen parameters. The calculations have been performed in case of 
transition metals Fe & Co for the wide range of compressions. It has been 
found that both the models yield close agreement with each other. The 
obtained results in the present paper satisfy the thermodynamic 
constraints in the limit of infinite pressure. 
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INTRODUCTION 

The Gruneisen parameter y is related to Thermoelastic 
properties of materials [1-3]. The pressure derivatives of y 
play a central role in predicting the thermal behavior , 
equation of state, and melting at high pressures [4-6]. The 
Gruneisen parameter y is an important physical quantity 
directly related to thermal and elastic properties of 
materials [7-9] as follows 
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Where a is the thermal expansivity K T and K s are 
isothermal and adiabatic bulk moduli, Cv and Cp are 
specific heats at constant volume and constant pressure, 
respectively. For a given material y decreases with 
increasing pressure P or decreasing volume V. 


On integrating Eq. (2), we get 
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And further integration gives 
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It has been found [5,8] that if we assume 
(A 0 — Aoo) = q 0 

Eq. (13) On differentiation yields 
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There have been various attempts [7-11] to formulate 
expressions for y(V) as well as y(P). The most critical 
examination of these formulations can be made by 
applying the thermodynamic constraints for the higher- 
order Gruneisen parameters [7,9, 15-17] at extreme 
compression in the limit of infinite pressure. In the 
presented study, we have been determined the higher 
order Gruneisen parameter in case of transition metals Fe, 
Co for the wide range of compressions. 

Formulation fory(V): 

According to the Stacey- Davis model we can write 
A = A M + (A 0 -A M )f (2) 

Qo 


In the present study we use Eq. (2), (3), (4) and (6) to 
obtain y, q, X and for transition metals. 


Formulation fory (P) 

The generalized free -volume formula for y as a function 
of pressure P, bulk modulus and its pressure derivative K 
= dK / dP has been widely used in literature [18-21]. This 
formula is written as follows [8-10] 


2 6 3V 


(i-vn) 


(7) 


where/is the free -volume parameter. Different values of/ 
such as 0, 1, 2 and 2.35 have been taken by Slater [18], 
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Dugdale and MacDonald [19], Volume dependence of Where y is a function of x such that y remains positive 
Griineisen parameter 199 finite at x-> 0 


Vashchenko and Zubarev [20], and by Barton and Stacey 
[21]. The main difficulty associated with Eq. (7] is 
regarding the uncertainties in the value of / arising from 
various approximations involved in the derivation of the 
formula [9]. It is even more uncertain to use Eq. (7] for 
evaluating derivatives of y because of difficulties in 
determining df / dP and d2 f / dP2 . Such calculations are 
tedious at finite pressures. However, it should be 
indicating that Eq. (7] in the limit of infinite pressure gives 
expressions [9, 15] which are free from the free -volume 
parameter f. An alternative formulation for_(P] has been 
presented by Shankar et al. [17]. They proposed the 
following relationship 



y 0 J K 


( 8 ) 


Where is the minimum value of K' in the limit of infinite 

pressure. Eq. (8] satisfies the following identity [20, 3] 

© =£ < 9 > 


Eq. (8] yields the following expressions for the higher 
order derivatives of the Griineisen parameter [16] 
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It should be mentioned that the expressions for y, q, X and 
given above satisfy the boundary conditions in the limit 
of infinite pressure viz. y^ and X^ both remain positive 
finite whereas qoo and xoo both tend to zero. Equations (8], 
(10], (11] and (12] at infinite pressure become identical to 
the corresponding expressions derived from the basic 
principles of calculus which have the status of identities 
[13,14,16]. 


INFINITE PRESSURE BEHAVIOUR 

P is a function of V such that increases continuously with 
the decreasing V. At extreme compression (V -> 0), P 
becomes infinitely large, then we can write 


dlnP\ 
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Since 

dlnP _ K 
dlnV ~ P 


Where K is bulk modulus 

, dP 
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Eq. (13] and (14] yield 



= positive finite 


Now taking y = - and x — V 


(13) 


(14) 

(15) 

(16) 


in expression 



= 0 


if y becomes zero at x-> 0 
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Because K/P is positive finite (Eq. 16], on solving Eq. (17], 
we get 

(1 - K P/ZCco = 0 (18) 


Eq. (18) yields 

= ©„ < 19 > 


Equation (19] has the status of an identity [7,5, 20]. In 
view of Eq. (16], it is revealed from Eq. (19] thatK the 
pressure derivative of bulk modulus, dK/dP, at infinite 
pressure remains positive finite. 


Result and discussion 

The quantities calculated as a function of pressure can be 
converted to those as a function of volume, and vice versa 
with the help of an equation of state theory. We use the 
Stacey K- primed equation of state [9] which is based on 
Eq.(9]. And can be expressed as 



On integration, Eq. (20] gives 
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Further integration of Eq. (21) gives 

O^-dO + t'" (H> 

The Griineisen parameter values y and its higher order 
derivatives q X and <f calculated from the Shanker 
formulation [16] using Eq. (8], (10], (11] and (12]. Values 
of y q X and get decrease with the increment in value of 
pressure. In the Stacey-Davis formulation they are 
expressed as a function of V/V 0 . At extreme compression, 
V tends to zero, q^ and both tends to zero, but 
7oo and X both remains positive finite. 


Conclusion 

The Griineisen parameter at extreme compression, the 
pressure for the earth inner core and outer core boundary, 
the melting temperature for Fe and Co is very close. The 
theory of materials at (oo) infinite pressure originally due 
to Stacey [5,7] and further developed by Shanker et al 
[12,13] has been extended in the present study. Values of 
y q X and <f at finite pressures for transition metals (Fe, 
Co] and Calculated from the two formulations compare 
with each other. Results are quite similar for different 
transition metals (Fe, Co]. Results can be produce using 
the input parameter given in table first. It should be 
emphasized that the pressure derivative of bulk modulus 
play the central role in determining Thermoelastic 
properties and equation of state of materials [1,3,21-23]. 
We thus obtain the following generalized expression for 
the higher pressure derivative of bulk modulus at infinite 
pressure. 
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Table (1) Values of input parameters for NaCl 



Z 

Yo 

Ko 

K 0 

K m 

wm 


A 0 


Fe(bcc) 

26 

1.83 

187.6 

5.8 

3.48 

1.33 

1.13 

6.13 

4.39 

Co(fcc) 

27 

1.77 

216.5 

4.7 

2.82 

1.24 

1.21 

5.91 

4.20 



- »T=300K 


- T=450K 


- T=600K 


Plots between pressure P and compression V/V 0 for Fe & Co at different temperatures 
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